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摘要 : 随 着 LAMOST 光谱 巡天 项 目的 不 断 推进 ,国际 上 积累 的 光谱 数量 已 经 步 入 干 万 量 级 。 
因 科学 研究 需求 或 大 量 低 信 品 比 光谱 的 不 断 产 出 , 在 光谱 处 理 过 程 中 , 常 涉 及 到 大 量 人 工 识 
谱 的 任务 。 为 方便 用 户 高 效 便捷 地 进行 人 工 识 谱 ， 本 文 基于 MATLAB 编程 技术 设计 并 实现 
了 一 套 一 维 光 谱 可 视 化 与 分 析 工 具 。 该 工具 包括 谱 线 标识 、 滤 波 去 噪 、 红 移 测量 、 等 值 宽度 
估算 等 功能 。 该 工具 操作 简单 、 直 观 ， 在 特殊 天 体 搜寻 ，LAMOST“UNKONWN ”光谱 的 
人 工 处 理 等 方面 具有 较 高 的 使 用 价值 。 
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Abstract: With the development of LAMOST spectroscopic Sky Survey project, the accumulated 
Spectra in the world have stepped into tens of millions of magnitudes. Due to the need of scientific 
research or the continuous output of a large number of low signal-to-noise ratio Spectra, a large 
number of manual Spectral identification tasks are often involved in Spectral processing. In order 
to facilitate users to carry out manual Spectrum recognition efficiently and conveniently, this paper 
designs and implements a set of one-dimensional Spectral visualization and analysis tools based on 
MATLAB programming technology. The tool includes Spectral line identification, filtering and 
denoising, redshift measurement, and equivalent width estimation. The tool is simple, intuitive 
operation. It has high application value in the search of special celestial objects and the artificial 
processing of LAMOST "UNKONWN" spectra. 
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1 引言 


随 着 各 个 大 型 光谱 巡天 项 目 不 断 地 开展 , 最 近 十 儿 年 来 , 天 文学 领域 已 收集 了 丰富 的 天 
体 光 谱 , 这 些 光谱 为 科学 研究 提供 了 重要 的 数据 基础 。 我 国 自主 研发 的 大 天 区 面积 多 目标 光 
纤 光 谱 天 文 望远镜 (Large Sky Area Multi-Object Fiber Spectroscopy Telescope，LAMOST)D 自 
2012 年 正式 巡天 到 2019 年 3 月 ， 已 经 获取 了 1125 余 万 条 光谱 ， 其 中 “UNKNOWN ”光谱 
100 余 万 条 。 随 着 巡天 观测 的 深入 ,“UNKONWN” 光 谱 会 越 来 越 多 。 
这 些 “UNKONWN” 光 谱 大 多 数 因为 信 吕 比较 低 ， 传 统 的 算法 无 法 处 理 ， 也 可 能 有 极 
少数 奇异 的 天 体 ， 由 于 谱 型 太 奇怪 ， 传 统 的 算法 无 法 识别 ， 如 毛 晓 燕 中 提出 了 基于 加 权 滤 波 
的 低 信 噪 比 LAMOST 光纤 光谱 信号 降 噪 。 由 于 每 条 光谱 的 获取 成 本 较 高 ， 且 这 些 光谱 本 身 
可 能 含有 丰富 的 科学 价值 ， 不 能 舍弃 它们 ， 因 此 如 何 识别 挖掘 这 些 “UNKONWN ”光谱 的 
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科学 价值 是 我 们 急需 解决 的 问题 ,一 些 天 文学 家 可 以 对 这 些 低 信 了 品 比 光谱 进行 分 类 和 参数 值 
估计 , 但 是 如 此 巨 量 的 数据 仅 靠 少 数 的 专家 利用 有 限 的 时 间 和 精力 完成 人 工 处 理 是 不 可 能 的 ， 
必须 开发 自动 化 的 光谱 分 析 软 件 。“ 工 欲 善 其 事 ， 必 先 利 其 器 ” 只 有 借助 良好 的 光谱 分 析 软 
件 ， 天 文学 家 才能 便捷 地 对 光谱 进行 分 析 和 研究 ， 从 而 提高 工作 效率 和 科学 产 出 。 

国内 外 专家 设计 了 诸多 一 维 光 谱 数 据 可 视 化 与 分 析 软 件 。 比 如 , 欧洲 空间 局 虚拟 天 文 台 
研发 团队 开发 了 VOSpecD， 可 以 对 光谱 进行 查询 、 分 析 、 拟 合 光 谱 能 量 分 布 ， 美 国 空间 望 
远 镜 科学 研究 所 开发 的 SpecView 钙 ， 可 以 对 光谱 进行 查询 、 可 视 化 与 分 析 ， 德国 虚拟 天 文 
台 团 队 与 捷克 共和 国 科学 院 的 天 文 研 究 所 共同 开发 的 SPLAT51， 可 以 显示 、 比 较 、 修 改 和 
分 析 天 文 光谱 数据 ， 并 且 能 查询 和 下 载 光 谱 ; 法国 天 体 物 理学 和 行星 学 研究 开发 的 CASSIS 
I， 可 以 查询 光谱 、 估 算 谱 线 的 柱 密 度 等 ， 中 国 科 学 院 国 家 天 文 台 开发 的 ASERA ， 可 以 
识别 类 星体 光谱 和 测量 红 移 ， 雷 国 洪 等 人 设计 了 一 套 基 于 数据 挖掘 的 专家 识 谱 平台 名 ， 可 以 
对 光谱 数据 进行 管理 、 可 视 化 与 分 析 和 科学 研究 与 科学 普及 等 。 这 些 软件 确实 能 挖掘 天 体 光 
谱 的 科学 价值 ， 然 而 我 们 发 现 VOSpec、SpecView、SPLAT 对 LAMOST 光谱 数据 存在 不 兼 
容 现 象 ， ASERA 主要 针对 类 星体 的 光谱 而 设计 ， 且 不 能 进行 滤波 去 噪 、 等 值 宽度 估算 等 。 

为 此 ， 本 文 基于 MATLABM 设 计 并 实现 一 套 一 维 光谱 的 可 视 化 与 分 析 工 具 ， 重 点 完成 
LAMOST“UNKNOWN” 的 人 工 处 理 ， 其 中 包括 光谱 的 人 工分 类 、 光 谱 的 红 移 和 等 值 宽度 
等 特征 物理 参量 的 测量 、 特 殊 天 体 搜寻 等 科学 目标 。 设计 的 可 视 化 界面 力求 简单 直观 , 适合 
普通 学 生 使 用 ， 通 过 发 动 学 生 力 量 ， 逐 渐 消 化 不 断 产 生 的 “UNKNOWN” 光 谱 数 据 。 
2 一 维 光谱 可 视 化 与 分 析 工 具 的 总 体 结构 设计 

本 文 将 设计 一 个 基于 MATLAB 开源 的 一 维 光谱 可 视 化 与 分 析 工 具 ， 利 用 MATLAB 设 
计 的 图 形 用 户 界面 MATLAB GUI) 实现 一 维 光 谱 的 可 视 化 ， 并 结合 相关 算法 实现 对 光谱 
进行 平滑 去 品 、 红 移 测 量 、 等 值 宽度 估算 等 分 析 功 能 。 一 维 光谱 可 视 化 与 分 析 工 具 的 总 体 结 
构 设 计 如 图 1 所 示 : 


图 1 一 维 光 谱 可 视 化 与 分 析 工 具 总 体 结构 设计 图 


Fig. 1 Overall architecture design of 1-D Spectral visualization and analysis tool 


一 维 光谱 可 视 化 与 分 析 工 具 分 为 三 个 部 分 ,分 别 为 数据 导入 , 数据 可 视 化 和 数据 分 析 。 
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(1) 数据 导入 : 用 户 将 本 地 光谱 数据 导入 到 一 维 光谱 可 视 化 与 分 析 工 具 中 ， 光 谱 文 件 
格式 主要 分 为 FITS 格式 、CSV 格式 和 ASCI 格式 。LAMOST，SDSS 等 大 型 巡天 项 目的 光 


谱 文 件 通常 为 FITS 格式 。 


标识 。 


(3) 数据 分 析 : 数据 分 析 包 括 红 移 或 视 向 速度 的 测 


估算 等 。 运用 定义 法 计算 红 移 或 视 向 速度 ， 运 月 


定义 法 测量 光谱 的 等 值 宽度 。 


3 一 维 光谱 可 视 化 与 分 析 工 具 的 具体 实现 


可 视 化 界面 设计 分 为 桌面 设置 、 光 谱 数 据 导入 设置 、 数 据 可 视 化 与 分 析 设 置 三 部 分 ， 数 
据 分 析 功 能 主要 包括 光谱 去 品 处 理 、 光 谱 红 移 或 视 向 速度 的 测量 和 等 值 宽度 测量 。 这 些 功能 


都 将 通过 Matlab 编程 技术 来 实现 。 
3.1 桌面 设置 
3.1.1 桌面 颜色 设置 


量 、 滤 波 降 品 处理 以 及 等 值 宽度 的 
上 小波 平 滑 去 品 算 法 进行 光谱 去 噪 处 理 ， 运 用 


(2) 数据 可 视 化 : 导入 光谱 数据 后 ， 通 过 MATLAB 程序 设计 实现 光谱 数据 可 视 化 。 
光谱 数据 可 视 化 包括 特征 谱 线 标识 、 多 光谱 加 县、 光谱 横 纵 坐标 轴 设 置 、 谱 线 类 型 设置 等 。 
标准 谱 线 文 件 是 在 静止 坐标 系 下 的 特征 谱 线 的 波长 列表 , 通过 导入 标准 谱 线 文件 来 进行 谱 线 


用 户 可 根据 个 人 的 喜好 选择 不 同 的 颜色 界面 ， 系 统 设计 了 6 种 颜色 可 供 选 择 ,分 别 为 绿 
色 (Summer Greens)、 蓝 色 〈Winter Blues)、 橙 色 〈Earthy Copper)、 粉 色 (Subtle Pinks )、 


灰色 (Gray Tones) 和 黑色 (Dark Blacks) 选项 。 


3.1.2 Logo 图 标 设 置 


3.2 数据 导入 设置 
光谱 数据 导入 过 程 如 图 2 所 示 : 


在 图 形 用 户 界 面 的 左上 角 是 系统 的 Logo 展示 区 ， 月 


国 select sor multiple fs spectra - 


日 户 可 设置 自己 的 Logo 网 标 。 
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ryfits 


DD spec-55892-F9202_sp01-143.fits 

DD spec-55920-F5592001_sp09-085.fits 
癌 spec-56597-GACO69N36B2_sp08-099... 
DD spec-57416-GAC100N13B1_sp06-046... 
DD spPLAT1489664908668692061.fits 

癌 spspec-51817-0399-235.fit 

DD spspec-52201-0692-378.fit 


~ [Spectra Fits Files (fits,~fid | 


oa (me) | 


图 2 FITS 光谱 文件 导 


入 过 程 


图 


Fig. 2 Diagram of importing FITS Spectral file 


系统 仅 文 持 导 入 一 维 光谱 数据 ， 
谱 数 据 上 传 至 系统 中 。 


3.3 数据 可 视 化 设置 
3.3.1 谱 线 标识 


UD 


主要 支持 FITS 文件 格式 ， 用 户 需要 将 本 地 FITS 光 


:201909.00191v1 


chinaXiIv 


ChinaXiv 合 作 期 刊 


谱 线 标识 能 辅助 用 户 识别 光谱 上 主要 的 特征 谱 线 , 在 静止 坐标 系 下 测 得 的 常见 的 特征 谱 
线 波长 如 表 1 所 示 : 


表 1 常见 特征 谱 线 的 波长 列表 (单位 :Nnm) 


Table 1 The wavelength list of some common Spectral lines (unit: Nnm) 


特殊 谱 线 | OvI Ly NV OI CI | Sily | CIY | HeI | OII 
静止 坐标 | 1033.8 | 1215.7 | 1240.8 | 1305.5 | 1335.3 | 1397.6 | 1549.5 | 1640.4 | 1665.8 
特殊 谱 线 | AT | CII | NelV | Mgll | NeV | NevI | Ol1l HeI SI 


静止 坐标 | 1857.4 | 1908.7 | 2439.5 | 2799.1 | 3346.8 | 3426.8 | 3727.1 | 3889 | 6718.3 


特殊 谱 线 | NI NI Li K H G MgI Nal CalT' 


静止 坐标 | 6529 | 6549.9 | 6707.9 | 3934.8 | 3969.6 | 4305.6 | 5176.7 | 5895.6 | 8500.4 
点 击 Speclines 按钮 ， 系 统 将 根据 用 户 设置 的 红 移 值 ， 将 特征 谱 线 以 虚线 的 形式 显示 出 


来 。 
3.3.2 多 光谱 加 赫 

点 击 Hold on 按钮 ， 图 形 用 户 界面 可 以 同时 导入 多 条 一 维 光谱 ， 并 对 这 些 光 谱 进 行 比较 
分 析 ， 点 击 clear 按钮 清除 可 视 化 界面 上 的 光谱 图 像 。 多 光谱 加 镭 显 示 如 图 3 所 示 : 


图 3 多 光谱 加 半 图 
Fig.3 Diagram of multi-spectrums overlapping 
3 中 可 视 化 界面 显示 了 两 条 不 同 的 类 星体 光谱 图 像 ,用 红色 虚线 显示 的 特征 谱 线 有 较 
强 的 吸收 线 , 特征 谱 线 MgII 发 生 红 移 后 波长 坐标 在 6500~6750nm 之 间 , 特征 谱 线 CIII 发 生 
红 移 后 波长 坐标 在 4500nm 左右 。 
数据 可 视 化 还 包括 对 光谱 图 像 进 行 网 格 划分 、 改变 光 谱 横 纵 坐 标 轴 的 范围 、 设置 谱 线 类 
型 和 设置 图 形 用 户 界面 的 标题 ， 用 户 可 以 根据 自己 的 需求 在 可 视 化 界面 的 编辑 框 进 行 设 


4 光谱 数据 分 析 功 能 的 实现 


根据 豆 广 民 no 等 人 的 小 波 去 噪 方法 ， 李 天 超 0 等 人 的 等 值 宽度 算法 以 及 刘 其 02 等 人 
的 红 移 测量 和 谱 线 证 认 方法 ， 本 文 运用 MATLAB 编程 技术 实现 了 光谱 平滑 去 噪 处 理 、 红 移 
测量 和 等 值 宽度 测量 等 光谱 分 析 功能 ， 
4.1 光谱 平滑 去 品 算 法 

小 波 去 品 方 法 既 能 去 除 图 像 中 的 噪音 ， 又 可 以 尽 可 能 地 保留 图 像 局 部 的 细节 ,因此 小 波 
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去 品 法 可 以 很 好 的 对 光谱 进行 平滑 去 噪 处 理 。 
4.1.1 光谱 平滑 去 噪 法 的 原理 
设 (X,y) 为 图 像 中 点 的 坐标 ，f (x,y) 为 含有 噪声 的 图 像 ，e(Xx, y) 为 网 像 上 噪音 干扰 ， 
o 是 噪声 方差 。 图 像 去 品 后 可 以 表示 为 : 
g(x,y)= f(x,y)—oe(x,y),(x=3700,3701,...,9000; y = —10,—9,...,10) 
小 波 去 噪 的 主要 步骤 : 
(1) 对 了 (x,y) 作 小 波 变换 ， 得 到 各 尺度 的 小 波 系数 wub ; 
(2) 设 定 各 分 解 层 阔 值 ， 对 小 波 系 数 W ;i) 进行 闵 值 处 理 ， 得 到 小 波 估计 系数 Wj ; 
(3) 重 构 小 波 系数 ， 再 利用 小 波 估计 系数 yt 有 进行 重新 构造 得 到 去 噪 处 理 后 的 图 像 


fly ys 
4.1.2 立 值 函数 选择 

阀 值 的 确定 和 对 小 波 系 数 的 阀 值 处 理 是 小 波 去 品 的 关键 , 这 里 选择 半 软 阀 值 法 进行 闪 值 
处 理 : 


Wjp) = 


、 sen (wo ) (wl-02) (Wop > 
0 (Wi < 14) 


其 中 ，wu 为 处 理 前 的 小 波 系数 ， 记 ;为 处 理 后 的 小 波 系数 ，? 为 阔 值 , 且 0<a <1。 


4.1.3 光谱 平滑 去 噪 结果 


本 文 将 光谱 平滑 程度 设置 成 6 个 等 级 ， 分 别 为 3、5、7、9、11 和 13〈 单 位 : pix)， 平 
滑 程度 为 3 pix 和 11 pix 的 光谱 图 去 噪 效果 如 图 4 所 示 : 
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图 4 光谱 去 噪 程度 比较 图 


Fig. 4 Comparisons of Spectral denoising 


从 上 图 可 以 看 出 , 平滑 程度 为 11pix 的 光谱 图 比 平滑 程度 为 3pix 的 光谱 图 更 平滑 , 视觉 
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效果 更 好 ， 但 是 11pix 可 能 会 损伤 部 分 微小 的 光谱 特征 ， 比 如 AIII 等 。 因 此 ， 用 户 需要 根 
据 自 己 的 需求 来 选择 平滑 程度 。 
4.2 红 移 或 视 向 速度 测量 算法 


在 物理 学 和 天 文学 领域 ， 红 移 现象 指 物体 的 电磁 辐射 由 于 某 种 原因 使 波长 被 拉 伸 的 现 
象 ， 具 体 表 现 为 光谱 的 谱 线 朝 红 端 移动 了 一 段 距 离 ， 即 波长 变 长 ， 频 率 降 低 。 红 移 现 象 目 前 


多 用 于 天 体 的 移动 以 及 规律 的 预测 上 ， 红 移 值 是 河 外 天 体 最 重要 的 物理 参数 之 一 。 
4.2.1 红 移 测量 法 的 原理 
设 谱 线 的 静止 波长 为 4 ， 观 测 波长 为 1，: 为 红 移 值 ， 则 
A4=(+z)A’ 
1 
= 二 -1 
Z 了 


谱 线 的 静止 波长 4 如 表 1 所 示 , 根据 可 视 化 界面 的 特殊 谱 线 位 置 确定 观测 波长 1 后 , 将 
= 观测 波长 输入 到 编辑 框 ， 点 击 enter 键 ， 得 到 特殊 谱 线 的 红 移 值 。 
ee 4.2.2 红 移 测量 结果 
本 文 以 三 条 光谱 测量 为 例 ， 运 用 一 维 光 谱 可 视 化 与 分 析 工 具 测 得 光谱 的 特征 谱 线 CIII 
和 MgI 的 红 移 值 如 表 2 所 示 : 


表 2 光谱 特征 谱 线 红 移 值 
Table 2 of Spectral redshift 


光谱 文件 名 CIII MgII 
spSpec-51817-0399-235.fit 1.3566 1.3675 
spSpec-52201-0692-378.fit 1.3570 1.3656 
spSpec-55920-F5592001-sp09-085 .fits 1.3587 1.3636 


4.3 等 值 宽度 测量 算法 


一 等 值 宽度 ， 即 与 吸收 《或 发 射 ) 谱 线 轮 廓 和 连续 谱 之 间 所 包围 的 面积 相当 的 高 度 为 1 
© 的 矩形 的 宽度 。 


4.3.1 等 值 宽度 测量 法 的 原理 
下 面 以 一 观测 光谱 上 CIII 谱 线 为 例 ， 概 述 天 体 特征 谱 线 等 值 宽度 的 测量 方法 (如 图 5 
所 示 )。 


图 5 CI 的 等 值 宽度 测量 示意 图 
Fig. 3 Chart of equivalent width measurement for CIII 


已 知 特征 谱 线 峰值 的 取 值 范围 是 [(wi, 九 ),(wL, 有 ,)] ， 峰 值 左边 点 的 取 值 范围 是 


(aa , 峰值 右 


峰值 右边 点 的 坐标 (wc2， 


边 点 的 取 值 范围 是 [(c,41),(ci,4i)] ,峰值 左边 点 


二 0 


wcl= 和 
J 
b, 

Pr : 
J 

fc2) 为 : 

We 
大 
d 

J 


(wl, fc; ) 是 坐标 (wcl, fel) 和 坐标 (wc2, fc2) 连 接线 上 的 一 点 ，(wl, fc; ) 的 值 为 : 


等 值 宽度 的 计算 结果 为 : 


4.3.2 等 值 宽度 测量 结果 


本 文 以 三 条 光谱 测量 为 例 ， 运 用 一 维 光 谱 可 视 化 与 分 析 工 上 其 测 得 光谱 的 特 包 


2— fcl 
fc; = Mw —wcl)+ fcl 


和 MgII 的 等 值 宽度 值 如 表 3 所 示 : 


表 3 光谱 特征 谱 线 等 值 宽 度 
Table 3 of Spectral Equivalent Width 


的 坐标 (ywcl fe1) 为 : 


E 谱 线 CII 


光谱 文件 名 CIII MgII 
spSpec-51817-0399-235 .fit 58.592 73.492 
spSpec-52201-0692-378 .fit 53.797 67.268 
spSpec-55920-F5592001-sp09-085 .fits 59.297 35.366 


5 总 结 与 展望 


LAMOST、2dF、SDSS 等 国际 上 大 型 巡天 项 目 不 断 获取 大 量 的 光谱 数据 ， 这 些 光 谱 为 


天 文学 研究 提供 了 重要 的 数据 基础 。 为 方便 
Matlab 研发 了 一 套 一 维 光 谱 可 视 化 与 分 析 工 具 。 该 工具 具备 以 下 几 个 方 


(1) 针对 一 维 光 谱 


谱 进 行 分 析 处 理 。 


(3) 能 像 SPLAT-VO 等 项 目 一 样 ， 实 现 


户 深 度 挖 掘 这 些 光 谱 的 科学 价值 ， 本 文 基于 
田 的 优势 


数据 ， 集 成 多 种 可 视 化 与 分 析 算 法 ， 方 便 用 户 简单 快捷 地 计算 出 光 


谱 的 红 移 或 视 向 速度 、 等 值 宽度 等 重要 物理 参数 。 
(2) 界面 简洁 ， 操 作 步 又 简单 ， 用 户 只 需要 掌握 可 视 化 界面 的 操作 方法 ， 就 可 以 对 光 


AS 


光谱 数据 可 视 化 与 分 析 ， 提 高 天 文学 家 光 
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谱 处 理 的 效率 。 
我 们 会 根据 需求 进一步 完善 该 工具 功能 , 借鉴 国际 上 优秀 的 光谱 可 视 化 软件 优势 , 最 终 
使 其 成 为 处 理 LAMOST、SDSS 等 光谱 数据 的 首选 工具 。 
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